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A Geometrical Method of Decoupling 

C. Baumgarten 1 

Paul Scherrer Institute, Switzerland^ 
(Dated: 7 September 2012) 

The computation of tunes and matched beam distributions are essential steps in the analysis of circular 
accelerators. If certain symmetries - like midplane symmetric - are present, then it is possible to treat 
the betatron motion in the horizontal, the vertical plane and (under certain circumstances) the longitudinal 
motion separately using the well-known Courant-Snyder theory, or to apply transformations that have been 
described previously as for instance the method of Teng and Edwards^. In a preceeding paper it has been 
shown that this method requires a modification for the treatment of isochronous cyclotrons with non-negligible 
space charge forces^. Unfortunately the modification was numerically not as stable as desired and it was still 
unclear, if the extension would work for all thinkable cases. Hence a systematic derivation of a more general 
treatment seemed advisable. 

In a second paper the author suggested the use of real Dirac matrices as basic tools to coupled linear optics 
and gave a straightforward recipe to decouple positive definite Hamiltonians with imaginary eigenvalues^. 
In this article this method is generalized and simplified in order to formulate a straightforward method to 
decouple Hamiltonian matrices with eigenvalues on the real and the imaginary axis. The decoupling of 
symplectic matrices which are exponentials of such Hamiltonian matrices can be deduced from this in a 
few steps. It is shown that this algebraic decoupling is closely related to a geometric "decoupling" by the 
orthogonalization of the vectors E, B and P, that were introduced with the so-called "electromechanical 
equivalence" . 

A mathematical analysis of the problem can be traced down to the task of a structure-preserving block- 
diagonalization of symplectic or Hamiltonian matrices, respectively. Structure preservation means in this 
context that the (sequence of) transformations must be symplectic and hence canonical. 

When used iteratively, the decoupling algorithm can also be applied to n-dimensional systems and is of 
order 0(n 2 ). 

PACS numbers: 45.20.Jj, 05.45.Xt, 41.85.-p, 03.30. +p 

Keywords: Hamiltonian mechanics, coupled oscillators, beam optics, Lorentz transformation 



I. INTRODUCTION 

The significance of the symplectic groups in Hamil- 
tonian dynamics has been emphasized for instance by 
A. Dragt 5 , and the Dirac matrices are known as gen- 
erators of the symplectic group 5p(4,i?) for long&. In 
Ref4 the author presented a toolbox for the treatment 
of two coupled harmonic oscillators that is based on the 
use of the real Dirac matrices (RDMs) as generators of 
the symplectic group Sp(4, R) and a systematic survey 
of symplectic transformations in two dimensions. This 
toolbox enabled to develope a straightforward recipe for 
the decoupling of positive definite two-dimensional har- 
monic oscillators. Here we present an improvement of 
the method that is based on geometric arguments, i.e. 
on the orthogonalization of 3-dimensional vectors asso- 
ciated via the electromechanical equivalence (EMEQ) to 
certain linear combinations of matrix elements. 

There is a long history of publications covering the 
diagonalization (and related) problems in linear alge- 
bra as well as in linear coupled optics, linear Hamilto- 
nian dynamics and control theory. A (non-exhaustive) 



list is given in the bibliography. Some of these publi- 
cations claim to have solved the problem of structure- 
preserving (block-) diagonalization and/or computation 
of the eigenvalues (see RefsJ£~— , but also Ref^ and 
references therein). However, none of the previous works 
(known to the author) takes full advantage of the group 
structure of the generators of 5p(4). The conceptually 
closest approach uses "quaternions" , the representations 
of which seems to be identical to the RDMs^i, but seems 
to be limited to orthogonal symplectic transformations. 
The decoupling method of Teng and Edwards has been 
the starting point, but turned out to be not general 
enough (see Refi^ and [D|) . 

The method that we present here, is based on a survey 
of all symplectic similarity transformations. We do not 
make specific assumptions about the Hamiltonian other 
than that it is a symmetric quadratic form and we present 
a geometric interpretation via the EMEQ, which provides 
a physical notation of otherwise complicated and mean- 
ingless algebraic expressions Q- Furthermore we believe 
that the presented Ansatz is an interesting example on 
how elements of classical physics, quantum mechanics, 
special relativity, electrodynamics, group theory, geomet- 
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ric algebra, statistics^ and last but not least symplectic 
theory fit together and allow to use a common formalism. 

The simplest classical linear dynamical system with 
interaction (coupling) has two degrees of freedom and 
hence a 4-dimensional phase space. It can be consid- 
ered as most "fundamental" and a detailled analysis of 
its properties will likely be instructive also for n > 2. 
Indeed it turns out, that the decoupling technique of the 
two-dimensional system can be iteratively applied to sys- 
tems with n degrees of freedom. A Jacobi-like iteration 
with pivot-search and is sketched in Sec. IV CI 

From the viewpoint of coupled linear optics the prob- 
lem is solved, if a symplectic transformation is derived 
that transforms (constant) Hamiltonian matrices to 2 x 2- 
block-diagonal form (see below). It has been shown in 
Ref.-, that the same transformation method can be ap- 
plied to symplectic matrices as well. The arguments will 
be shortly repeated below. When applied to symplectic 
matrices, the method is equivalent to the computation of 
the matrix logarithm. A solution for the couterpart, i.e. 
the computation of the matrix exponential with empha- 
sis on the use of Dirac matrices, has been presented by 
Barut, Zeni and Laufer in 19941. 

If 2 x 2-block-diagonal form has been achieved, the re- 
maining task is completely analogous to the application 
of the Courant-Snyder theory for one degree of freedom. 
Nevertheless some arguments require awareness of the 
eigenvalues and their relation to the properties of the 
Dirac matrices so that a reference to a complete diago- 
nalization seemed appropriate. 



II. COUPLED LINEAR OPTICS 

The Hamiltonian of a n-dimensional harmonic oscilla- 
tor with arbitrary coupling terms can be written in the 
form 



(1) 



where A is a symmetric matrix and "0 is a state-vector 
or "spinor" of the form V = (<?i,Pi, 92, P2, ■ • ■ , q n ,Pn) T - 
Even though the matrix A is time-dependent in the gen- 
eral case, it is well-known practice to use the Floquet- 
transformation to reduce it to a constant matrix for the 
treatment of periodic systems (see app. [B] and for in- 
stance Ref . 30 i 33 ). The symplectic unit matrix (usually 
labeled J or S) is an antisymmetric matrix that squares 
to the negative unit matrix. For n — 2 it is identified 
with the real Dirac matrix 70 . As described in Ref. 4 , 
it is possible to freely choose the order of the variables 
in the state vector. However, the order of the variables 
fixes the form of the symplectic unit matrix 70 0. We 



2 For n = 2, this also follows from the fundamental theorem of 
the Dirac matrices (see for instance Ref-JaL). 



prefer the use of a ordering system in which the phase 
space coordinates (qi,Pi) are grouped as pairs of canoni- 
cal conjugate variables, so that 70 has the form 



7o 



/ 1 

-1 





V 







1 

-1 0/ 



(2) 



Using the over-dot to indicate the derivative with re- 
spect to time (or pathlength), the equations of motion 
(EQOM) have the familiar form 



OH 



Pi 



dH 

'8q ( 



or in vector notation: 



^-7oV v ,H = Ft/. 
where the force matrix F is given as 



(3) 



(4) 



F = 7 oA. (5) 
From the definition of F one quickly finds t ha t 30 1 33 



F 1 = 70 F 70 . 



(6) 



where the superscript "T" denotes the transposed ma- 
trix. Matrices that obey Eqn. ([6]) are usually called 
"infinitesimally symplectic" or "Hamiltonian" 30 . Both 
terms are - in the opinion of the author - misleading: The 
former because F is neither symplectic nor is it infinites- 
imal, the latter since F does not appear in the Hamil- 
tonian while the symmetric matrix A does. In addition 
A and not F is in the view of the author the classical 
counterpart of the Hamiltonian operator (see App. A in 
Refi^). Furthermore, Eqn. © is a purely formal property 
and not necessarily connected to a Hamiltonian. There- 
fore the author decided to use the term "symplex" (plural 
"symplices") when referring to its formal definition (i.e. 
Eqn. ©) and its relation to the symplectic transfer ma- 
trix and to call it "force matrix" when referring to its 
physical content - especially with respect to the EMEQ 
(see Ref.— and below) . Accordingly we speak of an anti- 
symplex or cosymplex (i.e. "skew-Hamiltonian" matrix), 
if a matrix C fulfills the equation 



-70 C 70 . 



(7) 



If we write S (C) for (co-)symplices, respectively, option- 
ally with a subscript, then it is easy to prove that 



and 



Si S2 — S2 Si 

cs+sc, 



Sx S2 ~f~ S2 Sx 
Cx C2 ~h C2 Cx 

cs-sc 



symplex , 



(8) 



cosymplex . 



(9) 
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A. Dirac Matrices 

In the following we focuse on two degrees of freedom 
(n = 2), i.e. to a four-dimensional phase space and the 
use of the real Dirac matrices to describe its dynamics 
and transformation properties. Often the term "Dirac 
matrices" is used more restrictive and designates only 
four matrices, namely , k € [0 ... 3]. Here we consider 
the four basic Dirac matrices as the four basic elements 
of a Clifford algebra C7(3, 1) with 16 elements derived 
from the basic matrices (see app.[3}. For further details 
see for instance Ref.— i 31 i 32 . 

Any real 4 x 4-matrix M can be written as a linear 
combination of the RDMs 



15 



(10) 



fc=0 



The RDM-coefficients are given by0 



where Tr(X) is the trace of the matrix X. Only the first 
ten RDMs are symplices and since symplices obey the 
superposition principl o 4 ! 5 ' 33 , any force matrix (symplex) 
can be written as 



9 
fc=0 



(12) 



The solution of Eqn. ^ is known to be 

^( S )=exp(F S )V(0), (13) 
where the matrix 



M = exp (F s) 



(14) 



is called transfer matrix, which can be shown to fulfill 
the symplectic condition, if F is a symple x 4 i 5 i 33 : 



M 7o M = 7q . 



(15) 



Vice versa it is known that symplectic matrices can be 
written in the form of Eqn. ([T4P&2 3 -. 

Transfer matrices can be splittcd into two parts, one 
(M s ) being a symplex, the other (M c ) being a cosym- 
plesi^: 



M, 



(M- 7o M T 7o )/2 
(M + 7o M T 7o)/2, 



(16) 



which is in case of a symplectic matrix M identical to 

(17) 



M c = (M + M- 1 )/2 
M s = (M-M" 1 )^. 



3 Eqn. Hill is based on the fact that all RDMs except the unit 
matrix have zero trace. 



It has been shown in Ref4, that the decoupling of the 
symplex-part M s of a symplectic matrix M automati- 
cally decouples the corresponding cosymplex M c . Hence 
it is sufficient to derive a method to decouple symplices 
of the above mentioned type. In cases, where only the 
one-turn-transfer matrix is available, Eqn. (|16j) is used 
beforehand to extract the symplex-part of the transfer 
matrix. The decoupling algorithm can then be applied 
to this matrix (see also the detailed discussion in Ref4). 



III. BLOCK-DIAGONALIZATIOIM AND EIGENVALUES 

The force matrix F is by definition a product of a 
symmetric matrix A and of an antisymmetric matrix 70. 
Hence it has zero trace and the sum of all eigenvalues is 
zero. We restrict ourselves to systems with real- valued 
force matrices and therefore real-valued transfer matri- 
ces. The eigenvalues of real- valued 2 x 2-symplices are ei- 
ther both real or both purely imaginary (since the are the 
square root of a real expression). Block-diagonalization 
(in the case of the variable ordering as described above) 
means to find a symplectic similarity transformation R 
such that the matrix F = RFR 1 has the form 



Fi 
F 2 



(18) 



where F k are real 2 x 2-matrices. Since similarity trans- 
formations preserve the eigenvalues, a symplex is block- 
diagonalizable in the form that we are going to describe, 
if the (pairs of) eigenvalues are either real or imaginary. 
In case of imaginary eigenvalues, the corresponding de- 
gree of freedom (i.e. pair (qi,Pi)) is stable (or focused), 
while a pair of real eigenvalues belongs to an unstable 
(non-focused) degree of freedom. The corresponding be- 
tatron motion is unstable in the sense, that no sufficient 
focusing is present. 

However - in the general coupled case without further 
assumptions - F is a general 4 x 4-symplex (or larger). 
Using the RDMs it is relatively easy to construct matrices 
with complex eigenvalues. An example is 



F = E x 74 + B x 77 , 



(19) 



which has the complex eigenvalues ±i (B x ± i E x ). Since 
the eigenvalues are complex, also the 2 x 2-blocks are 
complex. They can be block-diagonalized, but the gen- 
eralization to the 2 n x 2 n-case requires a general treat- 
ment of the complex case, which goes beyond the scope 
of this paper. 

As in Ref 4- the author speaks of regular or massive sys- 
tems, if the Hamiltonian is positive definite and of irregu- 
lar or magnetic systems in case of indefinite Hamiltonian, 
respectively. Both types may be stable or instable and 
this distinction should not be confused with the question 
of stability. A detailed discussion of stability would go 
beyond the scope of this paper and we refer the reader 
for instance to Ref. 33 or Ref. ?7 and references therein. 
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A. The S-matrix 



IV. THE ELECTROMECHANICAL EQUIVALENCE 



The matrix of second moments a of a charged particle 
distribution 



has the time derivative 

cr = F a + a F 7 



(20) 



(21) 



Multiplication from the left with 70 and the use of 
Eqn. ([6]) leads to 



S = FS-SF : 

where the matrix S is defined by 
S = a 70 . 



(22) 



(23) 



If Eqn. ([23]) is compared to Eqn. ([5]), then it is obvi- 
ous that S is also a symplex as it is also the product 
of a symmetric and a skew-symmetric matrix and obeys 
Eq.© From Eqns. (H]), dHJ) and {20]) it follows that 



cr(s) = M(s)ct(0)M t (s) 



(24) 



The second moments of a matched distribution are un- 
changed after one turn (or sector) of period L so that 
cr(L) = er(0) so that one obtains in a few steps Q 



MS - SM = 0. 



B. The Eigensystems and Matching 



(25) 



Hence one finds that the matrices M, F and S have 
the same eigenvectors - but in general different eigenval- 



ue; 



,35.36. 



where 36 



F 

S 



A 
A 
D 



EAE" 1 
EDE 1 



M = EAE" 



(26) 



= Diag(e 



e — IU-l 



(27) 



= Diag(-iei,i£i, —ie2,ie 2 ) . 



oji are the oscillation frequencies and £i the emittances. 
If E is known, the second moments of the matched distri- 
bution can be computed by replacing the eigenfrequen- 
cies by the emittances. If a sympletic transformation R, 
is known, that brings F (and hence S and M) to block- 
diagonal form, then one can simply use the usual Courant 
Snyder theory for one-dimensional systems 3 -^. In this 
case an explicit computation of the eigenvectors is not 
required. 



It was shown in Ref.—, that the ten coefficients of the 
force matrix F or the S-matrix can be identified with 
energy £ and momentum P of & particle and with elec- 
tric and magnetic field (E and B, respectively) seen by a 
charged particle in external fields. The meaning of this 
identification is, that the corresponding coefficients of F 
or S transform under symplectic transformations in the 
exact same way as the fields and the momentum trans- 
form under the corresponding boosts and rotations. 

It was also shown that the envelope equations of cou- 
pled linear optics are isomorphic to the Lorentz force 
equation. The Lorentz group was found to be a subset 
of the two-dimensional symplectic group. The so defined 
"fields" (E and B) of the EMEQ should not be confused 
with the real fields of the beamline elements or accelera- 
tor components. 

This isomorphism has been named electromechani- 
cal equivalence (EMEQ). The ten possible symplectic 
transformations are identified with spatial and phase- 
rotations, Lorentz boosts and so-called "phase boosts". 
The transformation properties are analogous to those in 
Minkowski space-time. 

This structural analogy is the basic idea behind the 
electromechanical equivalence (EMEQ). Naturally, 70 is 
associated with the time-like components of 4-vectors 
(i.e. energy), the spatial matrices 7 = (71, 72, 73 ) T are 
associated with the momentum, the matrices 70 7 with 
the electric field and 714 70 7 with the magnetic field. The 
pseudoscalar has been named 714 = 70 71 72 73 (instead 
of 75, as convention in QED). The remaining six matri- 
ces are 710, which is the time-component of the pseudo- 
vector, (711, 712, 7i3) T = 7147 are the spatial compo- 
nents of the pseudo-vector and 715 = 1 is the unit ma- 
trix. A complete list is given in App. [A] further details 
in Ref A and in textbooks on quantum electrodynamics. 

The EMEQ is given by the following nomenclature: 



£ 
P 
E 
B 



fo 

(fi, h, fz) T 

(/•4, fb: h) T 
(JVj /8; /9) T : 



(28) 



with the f k given by Eqn. ([12]). Using the EMEQ, the 
eigenvalues of F (Eqn. [26] and Eqn. I2T[) can be expressed 
by: 



K 2 



UJ 2 

Det(F) 



£ 2 + B 2 - E 2 - P 2 

-2 £ P ■ (E x B) + £ 2 B 2 + E 2 P 2 

{E ■ P) 2 — (E ■ B) 2 — (P ■ B) 2 

{£ B + E x P) 2 - [E ■ B) 2 - (P ■ B) 2 



K 2 -4K 2 



4 See common textbooks on linear Hamiltonian dynamics or Ref™. 



(29) 
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Force matrices of stable systems have purely imaginary 
eigenvalues^!, so that for stable systems one has K 2 > 
and Kt > 2 y/K^. 

Using the notation of the EMEQ a general symplex F 
is explicitely given by 



F = 



E Z + By Ey — B Z 

E x —B x 



—E x 

E Z — By 

E y + B z B x E x 

-E, 



B x 

E Z — By 



— Ey — Bz 



—B x 

Pz £ P X 
S - P x Pz 



B z E z 




By 











Py 

-Pz £ + P x 
-£ + P x Pz 



(30) 



Note that F is block-diagonal, if B x = B z = E y = P y 
0, i.e. in vector notation: 



E ■ B 
PB 



(31) 





M r 


K 


Mi 


7o 


M r c + MgS 


M g c- M r s 


M b C 2 + E-^- S2 


7i 


M r C- {b) x S 


M g 


M b C- {r) x S 


72 


Mr C - (b)y S 


M g 


M b C-{r) y S 


73 


MrC- (b) z S 


M a 


M b C-{r)zS 


74 


Mr 


M g C+ {b) x S 


M b C+ (g) x S 


75 


Mr 


M g C+ {b)y S 


M b C+(g) y S 


76 


Mr 


M g C + (b) z S 


M b C+{g)zS 



TABLE I. Table of transformed "mass components" for 

symplectic transformations in 2 dimensions. Compare 
Eqns.[33[33]and|37] 



vector b equals the "Lorentz force" of a particle with 
magnetic charge, as the role of E and B is exchanged 
compared to g in the algebraic way that corresponds to 
a duality rotation about an angle of ~rA. 



One finds the following products: 



V. DECOUPLING OF 2-DIMENSIONAL SYSTEMS 
A. The geometrical approach 

In the following we describe a geometrical approach 
of decoupling that is inspired by the observation, that in 
the decoupled force matrix, the scalar products E ■ B and 
P ■ B vanish 4 . In Hamiltonian form (see Eqn. HIH below), 
also the product P ■ E is zero and only the components 
£ 1 P x , E z and B y remain. It is therefore instructive to 
analyze the symplectic transformation properties of these 
scalar products. The product E ■ B is known to be in- 
variant under rotations and Lorentz boosts. Formally it 
is a pseudo-scalar in contrast to the scalar component 
representing the mass. Hence one might loosely speak of 
"mass components" and use the abbreviations: 



M r = E-B 

Mg = B P 



M h 



E ■ P 



(32) 



The "mass components" are invariant under spatial ro- 
tations. We may therefore proceed with phase rotations 
and boosts. We introduce the following auxiliary vectors: 



r = EP + B x E 
g = £ E + P x B 
b = £B + Ex P, 

so that Ki from Eqn. [55] can be written as 



K 2 = P- M. 



Mi 



(33) 



(34) 



It is easy to see that g, f and b transform under spatial 
rotations just like usual vectors. It is also quite obvious 
that the vector g equals the usual Lorentz force and the 



9 = 



g ■ r 
g-b 
f ■ b 



-2£ P- {E x B) 
—2 £ P ■ (E x B) 
—2 £ P ■ (E x B) 
{£ 2 -B 2 )M b 
{£ 2 - P 2 ) M r 



£ 2 E 2 

£ 2p2 

£ 2 B 2 



B 2 P 2 
B 2 E 2 
E 2 P 2 



M 2 

9 

M 2 
M 2 



M r Mg 
MgM b 



(£ — Br) M g + M r M b 



(35) 

Furthermore we introduce the following abbreviations for 
a better readability 



K, 



and: 



Yi 
Zi 



c 

c 

C 2 



/•:, /•-', + n, Bj + 1\ Pj 

£ 2 - E 2 - B 2 + E 2 + Bf 
£ 2 -P 2 -B 2 + Ef + Bf 
£ 2 -P 2 -E 2 + Ef + Bf 



Pi 
Pf 
P? 



(36) 



cos (e) 
cos (2 e) 
cosh (e) 
cosh (2 e) 



s 
S 2 



sin (e) 
sin(2e) 
sinh (e) 
sinh (2 e) 



(37) 



The phase rotation generated by 70 yields: 

<f = gc + fs 
r 1 = fc — gs 
b' = b 



(38) 



The transformation of the mass components is listed in 
Tab. [H the transformed vectors g, f and b are given in 
Tab. [Til From the discussion of the normal form of the 
force matrix in Ref<^ it follows, that decoupling to block- 
diagonal form is done by a transformation that makes 
P y = E y = B x = B z = 0. Geometrically this means, that 
at least the two scalar products involving the magnetic 
field have to vanish. In consequence, the decoupling of 



6 



a two-dimensional harmonic oscillator requires the (par- 
tial) orthogonalization of the (3-dimensional) "vectors" 
E, B and P: 



M r 



E-B 



Mg = P B 




0, 



(39) 



which can be interpreted as a geometrical "decoupling" . 
A transformation to what we call "Hamiltonian" form 



/ a \ 
-P 
7 
\ -5 J 

requires additionally to make E x = P z = 0, i.e. 



(40) 



M h 



E-P^O. 



(41) 



The general form of symplectic transformations has been 
described in some detail in Ref^, here we give only a brief 
summary. A symplectic transformation matrix Rb is gen- 
erated by a symplex jb with b G [0 ... 9] and controlled 
by a parameter e: 



The last transformation is only possible, if \M r \ < \b y \ = 
\b[- 

(E-B) 2 < -2£P-(ExB) + S 2 B 2 +E 2 P 2 -(E-P) 

(44) 

The first transformations result in P ■ B = 0, so that 
Eq. 04] is identical to the requirement that > (see 
Eq. |29| . This means that the eigenvalues are either lo- 
cated on the real or imaginary axis, but not off-axis in 
the complex plane. If this condition is fulfilled, then the 
vector-components (g) y and (r) y , (b) x and (b) z are zero 
after the decoupling transformations have been applied. 

That is: we found a symplectic decoupling algorithm 
for both - systems with purely imaginary eigenvalues, 
which are called "strongly" stable^! and unfocused sys- 
tems with purely real eigenvalues. That the algorithm 
works in both cases equally well, is important for instance 
in the case of transverse-longitudinal coupling with space 
charge in cyclotrons^. 

We continue the discussion of force matrices with 
eigenvalues off axis in the complex plane in Sec. IV Bl and 
assume for now, that K2 > 0. Using the abbreviations 



Rb 



exp (7b §) 
exp(-7b§) 
Rb F R,, 1 



(42) 



M x = JM* + M* 



b% + b 2 



(45) 



The effect of a basic symplex 7^ depends on its "signa- 
ture" , which is positive for symmetric and negative for 
antisymmetric 7^: 

R f 1 cos (e/2) + 7 b sin (e/2) for 7fc 2 = -1 

\ 1 cosh (e/2) + 7fc sinh (e/2) for 7fc 2 = 1 ' 

where the bold printed 1 is the unity matrix. Note that 
transformations with = —1 (+1) are called rotations 
(boosts), respectively. The parameter e is called "angle" 
in case of rotations and "rapidity" in case of boosts. As 
the decoupling requires a sequence of transformations, we 
emphasize that the RDM-coefficients have to be updated 
according to Eq. [TT] after each transformation. 

Inspection of Tab.|I]shows that a straightforward strat- 
egy is the following: 

• M g — > 0: Make a phase rotation generated by 70 
with angle e = arctan(^-). This will always work 
independent on the size of Mj. 

• b — > \b\e y : Align the vector b along the y-axis 
by the spatial rotations with R 7 and an angle of 
e = arctan ( ) and with Rg about an angle of 

e = — arctan (I s -)- Such rotations can always be 
done. 



the RDM-coefficients of the block-diagonal (decoupled) 
force matrix are given by: 



p 



I _ P x Mr-Ey Mg Vb 2 -M2 



EL 



M r — > 0: Boost using 72 and angle , 



artanh( ; ^] 



b 2 (M r E x +M g P x )-£ b x Ml 
M x b yz \b\ 

jpl M r {by E z -b z Ey)+M g (by P * ~b * Py) 

Z ~ M x by z 

d/ _ £ B 2 -P-(ExB) 

y ~ \b\ 

B X = B Z = Ey = Py = 

(46) 

In order to bring the block-diagonal force matrix to 
Hamiltonian form, one may apply the following transfor- 
mations: 

• Mb — > 0: Use another phase rotation with 70 about 
e=i arctan (^) 

• P z — > 0: Use rotation about y-axis with 78 about 
e = — arctan (75^). 

After these two rotations, the matrix has Hamiltonian 
form, if K 2 > holds. In charged particle optics this is 
usually the case and therefore we consider this method 
as a generally applicable decoupling algorithm. 
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(r% 




7o 


(r% c - (g) x s 


(f) V c- (S)v s 


(r) z c- (g) z s 


7i 


(r) x C-M b S 


(>% 




72 


(r% 


(r) y C-M b S 




""M 




(rly 


(r) z C-M b S 


74 


(f% C 2 + % S2/2 


(r) y C + K xy S 


(r) z C + K xz S 


7« 


(7% C + K xy S 


{r) y C 2 + Y y S 2 /2 


(f) z C + K yz S 


76 


(r) x C + K xz S 


(r) v C + K yz S 


(r) z C2 + % S2/2 




13)'* 


i9)'v 


m 


To 


(g) x c+(r) x s 


(s)s c + (4 s 




7i 


(g) x C 2 ~ X x S2/2 


(g) y C-K xy S 


(g) z C-K XZ S 


72 


(g) x C-K xy S 


(g) y C 2 - Xy S2/2 


{g) z C-K zy S 


73 


{g) x C-K XZ S 


(g)y c - K yz S 


(g) z C 2 - X z S2/2 


74 


(g) x C + M b S 


(g) y 


(3). 


75 




(g) y C + M b S 


(g)z 


76 




(g) y 


(g) z C + M b S 




(py x 


(b)'y 


(b)' z 


7o 


(&)* 


(b)y 


(b) z 


7i 


{b) x C - M r S 


(b) y 


{b) z 


72 




(b) y C-MrS 


(b) z 


73 


(b) x 


(b)y 


(b) z C - Mr S 


74 


(b) x C + M g S 


(b)y 


(b) z 


75 




(b) y C + M g S 


{b) z 


76 




(b)y 


(b) z C + M g S 



TABLE II. Table of transformed vector components of g, f 
and b by symplectic transformations in 2 dimensions. Com- 
pare Eqns. 1551 1551 and 1571 



B. Complex Eigenvalues 

Even though the problem of complex eigenvalues has 
not yet been solved for the general 2 n x 2 n case, it is 
possible to give a solution for the 4x4 case as we are go- 
ing to describe here. The more general case of arbitrary 
2 n x 2 n-symplices with arbitrary (complex) eigenvalues 
can presumably be solved by a block-diagonalization with 
4 x 4-blocks for each set of complex conjugate eigenval- 
ues and 2 x 2-blocks for each pair of real or imaginary 
eigenvalues. 

If K2 < the eigenvalues of complex and a block- 
diagonalization with 2 x 2-blocks is not possible (within 
the reals). However a simplification of the matrix is possi- 
ble with the aim, that the RDM-coefficients of the trans- 
formed matrix have the following structure: 



P x = 

E x = 

E z ± 

£ = 







B x = 
E,, 



M n 



^ 
= 



Pz = 

B z = 

By ^ 

Mb = 0, 



(47) 



so that one finds g — and 6 = and the auxiliary vector 
r has only a single non- vanishing component r x . We 



distinguish two cases, the first with £ 2 < Max(P 2 ,E 2 ) 
and the second with £ 2 > Min(P 2 ,E 2 ). In both cases 
the goal is to let "energy" and "momentum" vanish by 
appropriate Lorentz or phase boosts. Then one may align 
B along the y-axis and rotate about the y-axis to make 
E x = 0. Then the conditions of Eqn. HTJare fulfilled. 



1. The Low Energy Case 

The decoupling strategy for the first case, i.e. for £ 2 < 
M&x(P 2 ,E 2 ): 

• M g — > 0: Apply a phase rotation Ro with angle 
e\ = arctan(^-). Note that this maximizes M r = 
E ■ B. 

• E — > I IS I e y : Align the vector E along the y-axis 
by the spatial rotations with R7 and an angle of 
£2 = arctan (-f^) and (after an update of the RDM- 
coefficients and a recomputation of the auxiliary 
vector and mass components) with Rg about an 
angle of £3 = — axctan (jf-)- 

• £ — > 0: Boost using R2 and rapidity £4 = 
artanh(-J-). According to the assumptions, this is 
possible and does not change E x — or E z — 0. 



• P T — > 0: Boost using R 3 and rapidity £5 = 



- art anh (-g^) 



• P z 0: Boost using Ri and rapidity Sq = 
artanh(-g^). Since £ = E z — E x = 0, the energy £ 

as well as E are unchanged by the boost. 



B — > \B \ e y : Align the vector B along the y-axis 
by the spatial rotations with R7 and an angle of 
£7 = arctan {^-) and (after an update of the RDM- 
coefficients and a recomputation of the auxiliary 
vector and mass components) with Rg about an 
angle of e$ = — arctan ( : § 2L ). 



• EL 



0: 



Rotate about y-axis with Rs with an 

E, 



angle of £g = arctan (§ t ). 



2. The Intermediate Energy Case 

The case where £ 2 > Min(P 2 ,£; 2 ) but K 2 < might 
be called "intermediate", since the energy is large com- 
pared to the "low energy" case, but not large enough to 
make K 2 > 0. The following procedure lead to the state 
described by Eqn. 27] 

• Mh — > 0: Apply a phase rotation Ro with angle 
£1 = i arctan ( Sr^h, )■ Note that this transforma- 

tion minimizes P 2 . 
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P — > \P\ e y : Align the vector P along the y-axis 
by the spatial rotations with R7 and an angle of 
£2 = arctan (^p 5 -) and (after an update of the RDM- 
coefficients and a recomputation of the auxiliary 
vector and mass components) with Rg about an 
angle of e 3 = — arctan (-p^). Since M& = E ■ P = 
one also has now E y = 0. 

P y — > 0: Lorentz boost using R5 and rapidity £4 = 
artanh(-y-). 

B —> \B \ e y : Align the vector B along the y-axis 
by the spatial rotations with R7 and an angle of 
£6 = arctan (J 2 -) and (after an update of the RDM- 
coefficients and a recomputation of the auxiliary 
vector and mass components) with Rg about an 



angle of £7 = — arctan 



)• 



£ ->• 0: Boost 
artanh(-J-). 



using R2 and rapidity £§ = 



• E x — > 0: Rotate about y-axis with Rs with an 
angle of £g = arctan (-f^). 

In both cases the transformed matrix F then has the 
form 



F = E y 75 + E z 7 6 + B y 7s . 



(48) 



In order to bring it to Hamiltonian form, one applies the 
transformation R2 with an "angle" of itt/2: 



R 



1 



exp(i7r/47 2 ) = -j= (l + «7 2 ) 



so that F RF R 1 is: 



-i E y 70 + E z 7 6 + B y 7s . 



(49) 



(50) 



Note that the complex eigenvalues of a force matrix with 
K 2 < all lie on a circle of radius p = (Kf + 4 \K 2 \) 1/4 
in the complex plane. 



C. Decoupling n-dimensional Symplices 

The general eigenvalue problem of symplices (Hamil- 
tonian matrices) is an area of intense research. The al- 
gorithm presented above is based on a physical and geo- 
metrical analysis of 2-dimensional linear symplectic sys- 
tems. As described before, the algorithm is limited to 
symplices that have real or imaginary eigenvalues, but a 
generalization to include complex eigenvalues might be 
possible - even though not urgently required in charged 
particle optics 0. 



5 The S-matrix for instance will never have complex eigenvalues 
as it is derived from the matrix of second moments. Complex 
eigenvalues would only be possible for correlations with a mod- 
ulus greater than one. 
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FIG. 1. Solid line: Number of iterations required to bring a 
2nx2n symplex (Hamiltonian matrix) to Hamiltonian form. 
Dashed line: Approximation by 5 " ^ ■ The number rib of 
non-diagonal 2 x 2-blocks is nt = - ^" 2 ~ 1 ' > ■ 



In order to decouple symplectic systems with more 
than two degrees of freedom, the described algorithm can 
be used in an iterative scheme analogous to the Jacobi 
method for symmetric matrices. If all eigenvalues are 
real or imaginary, it is possible to avoid computations 
using complex numbers. The 2 n x 2 rt-symplex is then 
regarded asanxn matrix of 2 x 2-blocks. A possible pivot 
search picks the maximum average square amplitudes of 
all non-diagonal blocks and starts the decoupling with 
those blocks. The iterative scheme allows to compute si- 
multaneously the symplectic transformation matrix and 
the resulting block-diagonal or Hamiltonian form sym- 
plex with high precision. 

Given {x} is a sequence of random numbers between 
zero and one, then one may construct random symmetric 
2 n x 2 n-matrices A according to the rule: 



x — h for i j 
n + x for i = j 



(51) 



The increase of the diagonal terms helps to avoid complex 
eigenvalues. The symplex to decouple is then given by 
F = 70 A. We tested the algorithm with these random 
matrices up to n = 12 and logged the number of 4 x 4- 
diagonalization steps. Fig.[T]shows the average number of 
iterations that is required to compute the transformation 
that brings a 2 n x 2 n symplex to Hamiltonian form, i.e. 
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into the form: 





-7i 





V o 



D. Diagonalization 










o \ 


o / 



(52) 



In order to proceed from Eqn. 00] towards diagonal- 
ization, the matrix is scaled using the generators 73 and 
74: 



R 

s 
t 



74)f] 



exp [(73 + 74)5 + (73 
Diag(exp (— s), exp (s), exp (— t), exp (t)) 
exp (log (|f |)/4) 

1*1 



(53) 



exp(log(|2|)/4) 



so that one obtains (for stable systems), what we call 
"normal" form: 



(54) 



where the signs of the frequencies u)\ and L02 can both be 
positive and negative, depending on the signs of a, (3, 7 
and 5. At this stage, all components of g and r as well 
as (b) x and (b) z are zero. Only (b) y is non-zero. 



RFR 


-1 






( ° 







\ 






















UJ 2 


V 





-LJ 2 


/ 



E 



E 7oE^ 

En 1 



Eo Diag(iwi, - 
( l-i -1 
l + i 



l-7o 



V 



§( 

70 



l + i 



«73 - 



iWl, !W 2 

H 


1 -i 
l + i 
»7e) 



-iw 2 )E 1 

\ 


-l + i 
l + i J 



E, 



That is - the last transformation matrix that is required 
for diagonalization is not only symplectic - it is also uni- 
tary. 



E. Example 

A simplified and idealized cyclotron model with space 
charge was described, which served as an example for an 
irregular system 3 ' 1 *. Without repeating all details, the 



constant force matrix has the following form: 

/ 1 \ 

-k x + K x h 
-h ^ 
V K zl 2 J 

The RDM-coefficients are then given by 
£ = 

Px = 

Py - 

E x = 

Ey = 

E z = 
B„ 



I'Voin this one finds for the 
g, r and b: 



(56) 




(57) 



'mass" terms and the vectors 



M r 

M b 

9 
f 
b 



E-B = -\(\ + K zl 2 ) 
P B = 
E- P = 

(o,f(if 27 2 -i),i±^^r 



(58) 



( 
(0 



4 7 2 



L ,0,0) J 

4 I 



According to the geometrical approach, the first trans- 
formation can be omitted, since the "mass" M g is zero. 
The second transformation using 77 aligns b along the 
y-axis. The second rotation may again be omitted, since 
the vector r is already aligned along the x-axis. The 
last transformation is a phase boost using 72 and is suffi- 
cient to bring F into block-diagonal form. This transfor- 
mation would usually change the value of Mb, but here 
it does not, since Mf, = (r) y = as can be seen from 
Tab. [U Hence M b remains zero - M g is invariant under 
both transformations. Hence, all "mass terms" are then 
zero after the described two transformations so that the 
system is decoupled. 



F. Operators, Expectation Values and Lax Pairs 

Coupled linear optics is in its essence (as quantum me- 
chanics) a statistical theory. Since the reference trajec- 
tory is fixed, the coordinates are always taken relative 
to the local reference frame and the geometry is (only) 
locally euclidean. Even though the starting point is the 
description of single particle motion, the orbits of sin- 
gle particles are usually both, hard to access experimen- 
tally and of low practical value. The description of the 
beam by average values in contrast is both - measureable 
and of high value. The use of symplectic transformations 
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leaves the expectation values unchanged. We can there- 
fore evaluate the expectation values of any operator O in 
an arbitrary reference frame: 



(O) = (^OV> 

= (Vj T 7 R _1 ROR- 1 R^) 

= (V T R T 7o6^) 

since for symplectic R we have 

R T 7o = 70 R" 1 
R T 7o R = 70 



(59) 



(60) 



The time derivative of the expectation value of an ar- 
bitrary operator O, that does not explicitely depend on 
time, is: 



= V' T F T 7 O OV' + ^OF-0 
= ^(OF-FO)tA 



Equations of the form (here S = o"7o) 
S = FS - SF, 



(61) 



(62) 



appear frequently in the theory of coupled linear optics 
and it is worth mentioning that Eq. [52] is a so-called Lax 
representation and the operators S and F are a so-called 
Lax paii^^. As a consequence, the expressions 



h = Tr(S k ) 



(63) 



are first integrals of motion, where TVQ is the trace. Us- 
ing again the EMEQ to express the elements of S, one 
finds: 

h = Tr(S) = 

h = Tr(S 2 ) = -4 [s 2 -P 2 +B 2 - £ 2 ) = -4 K x 

h = Tr(S 3 ) = 

h = Tr(S i )=A(K 2 +4K 2 ) 

(64) 

The values of K\ and K 2 are (as expected) first integrals 
and constants of motion. The complete expression for S 4 
is 

S 4 = (K 2 +AK 2 )1-AK X (M g 710 + M r 7l4 + b 714 7 

(65) 

Another derivation of Eqn. [53] has been given in 4 ^. 

VI. SUMMARY AND OUTLOOK 

A powerful method for symplectic decoupling of the n- 
dimensional non-dissipative harmonic oscillator has been 



developed. The method apparently is of order 0(n 2 ). 
The resulting block-diagonal symplex can be used to 
compute the cr-matrix of matched beam ellipsoids of lin- 
ear coupled systems in charged particle optic a 3 ' 4 ' 36 . The 
presented parametrization gives deep insight into the 
general nature of coupling and might be instructive also 
in other areas of physics. The algebraic problem of find- 
ing the eigenvalues and eigenvectors of a two-dimensional 
symplectic system was solved using geometrical argu- 
ments based on the use of the real Dirac matrizes and 
the electromechanical equivalence. 
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Appendix A: The 7-Matrices 

To complete the list of the real 7-matrices used 
throughout this paper: 



70 = 



72 

714 
74 
75 

76 

710 
711 
712 
713 



/ 1 \ 

-10 
1 

V -1 j 
I 1 \ 

10 
10 

V 1 J 
70 71 72 73; 
7071; 
7072; 
70 73; 
71470 
71471 
71472 
71473 



71 = 



73 



/ 



715 
77 

78 

79 

71 72 73 
70 72 73 
70 73 71 
70 71 72 



/ 

-1 



V 

-1 




V 
1 

714 70 71 

714 70 72 
714 70 73 



-1 N 

1 
10/ 
N 



1 
0-10 
1 / 

' 72 73 
: 73 71 
: 71 72 



(Al) 



Appendix B: Floquet Theorem 

If the matrix A in Eqn. (TTJ) and hence the forces are not 
constant, but periodic (F(t + T) = F(t)), then Floquet's 
theorem can be applied and the solution has the general 
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form^^: 

M(t) = K(t) exp(Fi), (Bl) 
where K(t) is symplectic and periodic with period T. 



M(0) = 1 => K(0) = 1 
K(t + T) = K(i) K(T) = 1 



(B2) 



The transfer matrix of one period of length T ( "one-turn- 
transfer-matrix" ) 



M(T) = M T =exp (FT) 



(B3) 



is identical to the transfer matrix for a system with the 
constant force matrix F and length T. In this sense F is 
the "average" or "effective" force matrix with respect to 
one turn and can formally be written as^: 

F = I In (M T ) . (B4) 

From Eqs. IB II one derives in a few step o 30 ^ 38 : 

K = F K K F . (B5) 

If the canonical transformation represented by K has 
been applied to the state vector, then with K(0) = 1 
it follows: 

i){t) = K - 1 ip(t) = exp (F f) -0(0) = M ip(0) . (B6) 

Note that the knowledge of K is not required to solve 
the decoupling problem, as long as the one-turn-transfer 
matrix Mt is known. Mt can either be obtained as a 
product of the transfer matrices of all beamline elements 
or simply by numerical integration. If the matched beam 
distribution has been found at an arbitrary (known) po- 
sition s = along the closed reference orbit, then the 
matched distribution can be computed for any position 
s using: 



a(s) = M(s)a(0)M T (s) 



Appendix C: Quick Guide to Decoupling 



(B7) 



To start with it is required to have either the average or 
constant force matrix F or the symplectic transfer matrix 
M that represents a complete turn or (cyclotron) sector. 
In the latter case one computes an auxiliary force matrix 
by 



M s = i(M + 7o M T 7o ), 



(CI) 



while the usual (effective) force matrix has the form 

F = EDiag(« uji, -iui,iu)2, —icJz) E -1 , (C2) 

uii being the betatron frequencies. The auxiliary matrix 
has the same structure 

M s = EDiag(isi,-isi,is 2 ,-is 2 )E" 1 , (C3) 



but different eigenvalues Sj = sin (uji r) , where r 
27r Qi with the betatron tunes Qi. 

Now compute the RDM-coefficients according to: 



£ 


= -Tr(F 7o 


+ 7o F)/8 


P, 


= Tr(F 7l + 


7iF)/8 


Py 


= Tr(F 72 + 


72 F)/8 


P, 


= Tr(F 73 + 


73 F)/8 


E T 


= Tr(F 74 + 


74 F)/8 


E y 


= Tr(F 75 + 


7 5 F)/8 


E z 


= Tr(F 76 + 


76 F)/8 


B x 


= -Tr(F 77 


+ 77 F)/8 


By 


= -Tr(F 78 


+ 78 F)/8 


B z 


= -Tr(F 79 


+ 79 F)/8 



(C4) 



Note that the coefficients for 7 fc with k £ [10, . . . , 15] 
must be zero - otherwise the system is not symplec- 
tic. Then compute the eigenvalues and auxiliary vectors 
f, g, b according to Eq. [29l [32] and [33] Construct the 
transformation matrices Rb according to: 



R 



b — 



1 cos e/2 + 7 b sine/2 for b G [0, 7, 8, 9] 
1 coshe/2 + 7b sinhe/2 for 6e[l,...,6] 

1 cose/2 - 76 sine/2 for b e [0, 7, 8, 9] 
1 cosh e/2 — 7 6 sinhe/2 for 6e[l,...,6] 

(C5) 

Transform with 70 and e = arctan {^m^ '■ 



F > Rq F R n 



(C6) 



Recompute RDM-coefficients, then transform using 77 
with e = arctan (jrj ■ Recompute RDM-coefficients, 

then transform using 79 with e = arctan^ ^- ^. Re- 
compute RDM-coefficients, then transform using 72 with 
e = artanh^-^^. The (auxiliary) force matrix should 
now be block-diagonal. Recompute RDM-coefficients, 
then transform with 70 and e = arctan ( g 2 M p 2 j . Re- 
compute RDM-coefficients, then transform with 7§ and 
e = — arctan ^j^ 5 -^ . Now the (auxiliary) force matrix 

should have normal form, so that the frequencies (or their 
sines) are given by: 



>^2 



The complete transformation is given by: 



(C7) 



Rn ' Ri 



■ R-n 1 

■ Ro 



(C8) 



R 1 

R = Rn • Rn-1 • 

F d = RFR 1 

If the auxiliary matrix has been used, then compute the 
matrix M c according to 

1 



Mr 



■R(M- 7o M T 7o)R- 1 



(C9) 
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The cosines of the tunes are then given by: 



Appendix E: Cosymplices 



cos (uii t) + cos (w 2 r) = Tr(M)/2 
cos (wi r) — cos (u>2 t) — Tr (M 712 



712 M)/4 

(CIO) 



Appendix D: The Teng and Edwards Ansatz 

Assume that we have an even number of DOF, so that 
the 4n x 4n transfer matrix R can be written in block- 
form according toil 2 -: 



R 



M m 
n N 



(Dl) 



where all quadratic submatrices are of size 2 n x 2 n, then 
the matrix R is symplectic, if 



70 




11170 m 
N7 m T 



M 70 n T + m 70 N T 
^ T 1 ■ N 70 N T 



n7on J 



which yields: 



To 
To 




M 7o M J 
1170 n T ^ 
M7 n T 



+ m 70 nr 
N 70 N T 
f m 7o N T . 



(D2) 



(D3) 



where 70 has - in dependence of the context - to be taken 
as2nx2nor4nx4rt. 

If one now assumes that M = N = C 1, then it follows 
that 



7b (1 
7b (1 



C 2 ) 
C 2 ) 



n = 



11170 m 
n7 n T 

70 m T 70 



(D4) 



= m . 



If one assumes with Teng and Edwards, that C = cos ((f)), 
then may define m = sin ((f)) m s and n = sin ((f)) n s with 
symplectic matrizes m s and n s , respectively: 



To = 

To = n s 7 n 



m s 70 

T 



(D5) 



It has been shown in Ref. 3 , that C — cos ((f)) is not the 
general case, since one might also choose C — cosh (4>), 
m = sinh (0) m s and n = sinh(</>)n s . In this case one 
finds 



-To 
-To 



m s 70 m; 

n s 70 , 



(D6) 



i.e. the matrizes m s and n s can also be antisymplec- 
tic (symplectic with multiplier —1). Still the matrix R 
remains symplectic. 



The geometric approach is based on the second order 
terms, i.e. products of the RDM coefficients. It is there- 
fore instructive to see where else these terms appear. For 
instance one quickly finds the "mass" terms and vectors 
g, r, b in the following products: 



FF 


= -(£ 2 -P 2 + B 2 -E 2 )l + 

+ 2 Mg 710 +2 6714 7 


2 M r 714 


F 7o F 


= {3£ 2 -P 2 -E 2 -B 2 ) l0 - 
+ 2r*7 + 2(7 7o7 + 2& 714 To T 


4£F 


F714F 


= 2 M b 710 -2M r l 






+ (£ 2 -P 2 + E 2 -B 2 ) 114 + 


2 3T14T 


F710F 


= 2 M b 714 - 2 M g 1 






+ (£ 2 + P 2 -E 2 -S 2 ) 7 io + 


2 r 7147 



(El) 



So that in the decoupled and normalized case (see 
Eq. [54]), these products are: 



-(£ 2 -P 2 + B 2 -E 2 )l + 2 (b) y 712 

(3£ 2 -P 2 -E 2 -B 2 ) l0 -A£~F 

2 (b) y 78 (E2) 

(£ 2 -P 2 + E 2 -B 2 ) lu 

(£ 2 + P 2 -E 2 -B 2 ) lw 



FF 

F 7o F 

F714F 
F710F 



Appendix F: Expectation Values (Complement) 

In Ref. 4 it has been shown that the expectation values 
of the RDMs, f k , defined by 



fk 



1 



■1p Ik ip , 



(Fl) 



vanish for all cosymplices, i.e. for 7^ with k € 
[10, . . . , 15]. It was also shown that for all symplices (i.e. 
7^ with k € [0, . . . , 9] or linear combinations thereof) 
the expectation values gt = ^(7fc F + F7fc)^ vanish. 
Nevertheless nothing was mentioned about the gk for 
k E [10, ...,15]. The complement is given in the fol- 
lowing: 



.910 = 2 (P x f 7 + P y f 8 -\ 

.911 = 2 (-£ f 7 + B x f Q 

.912 = 2 (-£ h + B y fo 

913 

9U 



2 (-£ h + B z f 
2 {E x fr + E y f s - 



.915 = 2(f; 



Pz h -B x h- By f 2 + B z h) 

- P,. h +Eyh~ Py h ~ E z f 2 ) 

-P x f 6 + E z h -P Z U- E x h) 

-Pyfi + E x f 2 -P x f 5 -E y f 1 ) 

- E z fg — B x fi — By f 5 + B z fo) 

(F2) 



According to Eq. [in] the expectation values of the opera- 
tors gk are: 



g k =^( 7 fcF 2 -F 2 7fc )V. 



(F3) 



13 



The square of the force matrix is given in Eq. lEll Now we 
insert this into Eq. IF3I The scalar part commutes with 
all 7^ and hence contributes nothing. Since all commuta- 
tors of symplices with cosymplices result cosymplices, we 
obtain gk = for all symplices. This had to be expected 
as for all symplices we had gk — 0. Hence the remaining 
terms are: 





= 4^(M r70 4 










= 4(M r / + 6 


x J4 + by J 5 4 


-b z 


h) 


5n 


- 2(41/,./!- 


AM g U + by 


h 
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